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ABSTRACT 


A program capable of calculating the deBign-point and quasi-steady-state 
annual performance of a paraboloidal-concentrator solar thermal power plant 
without energy storage has been written for a programmable calculator equipped 
with suitable printer. The power plant may be located at any site for which a 
histogram of annual direct normal insolation is available. 

Inputs required by the program are aperture area and the design and 
annual efficiencies of the concentrator; the intercept factor and apparent 
absorptance of the receiver aperture, and the receiver heat loss; the design 
efficiency of the power conversion subsystem and a polynomial representation 
of its normalized part-load efficiency; the efficiency of the electrical 
generator or alternator; the efficiency of the electric power conditioning and 
transport subsystem; and tl e fractional parasitic losses for the plant. 

(Losses to auxiliaries associated with each individual module are to be 
deducted when the power conversion subsystem efficiencies are calculated.) 

Outputs provided by the program are the system design efficiency, the 
annualized receiver efficiency, the annualized power conversion subsystem 
efficiency, the total annual direct normal insolation received per unit area 
of concentrator aperture, and the system annual efficiency. 
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SECTION I 


INTRODUCTION 


A program capable of calculating the design-point and quasi-steady-state 
annual performance of a paraboloidal-concentrator solar thermal power plant 
operating without energy storage is described in this report. The program has 
been written for a programmable calculator equipped with a suitable printer, 
viz., the Texas Instruments TI-59 with PC-100C. For computational purposes, 
the power plant may be located at any site for which a suitable histogram of 
annual direct normal insolation is available. 

Inputs required by the program are the design and annual efficiencies of 
the concentrator; the intercept factor and apparent absorptance of the 
receiver aperture, and the receiver heat loss per unit area of concentrator 
aperture; the design efficiency of the power conversion subsystem and a 
polynomial representation of its normalized part-load efficiency; the 
efficiency of the electrical generator or alternator; the efficiency of the 
electric power conditioning and transport subsystem; and the fractional 
parasitic losses for the plant. (Losses to auxiliaries associated with each 
individual module are to be deducted when the power conversion subsystem 
efficiencies are calculated.) 

Outputs provided by the program are the system design efficiency, the 
annualized receiver efficiency, the annualized power conversion subsystem 
efficiency, the total annual direct normal insolation received per unit area 
of concentrator aperture, and the system annual efficiency. In addition, 
pertinent performance data for the concentrator, receiver, and power 
conversion subsystems, and for the entire system are output at each of the 
twenty-one median values of direct normal insolation included in the histogram 
(0.025, 0.075, . . ., 0.975, 1.025 kW/m 2 ). 

Input parameters required for the calculation of paraboloidal-dish solar 
thermal power plant performance, together with supporting explanatory notes, 
can be entered into a form such as the one presented in Table 1. This form 
also provides spaces for listing all the design and annual performance data 
pertinent to a given configuration operating under given design and annual 
insolation conditions. 

The input and output variables and a description of the program and 
corresponding equations are presented in Section II, Description of the 
Problem and Its Solution. A flow chart and thoroughly annotated program 
listing are included in Appendixes A and B, respectively. 

A sample problem statement and corresponding solution are presented in 
Section III, Sample Problem and Solution; input data for the sample 
problem are included in Appendix C. 

Section IV, entitled User Instructions, provides step-by-step 
instructions for using the program. 

Section V provides the concluding remarks. 
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Table 1. Performance Summary Input Data Form 


Technology: 

Time Frame: 

2 

Design Direct Normal Insolation: kW/m 

Geographic Location: 

Annual Direct Normal Insolation: kWh/m^/y 

Concentrator 
Design Parameters 

Aperture area (m 2 ) 

Reflectivity . 

Blocking . . 

Efficiency . 

Annual Operating Parameters 

Degradation 

Shading 

Efficiency 

Receiver 

Design Parameters 

Aperture area (m 2 ) . . . 

Temperature of the thermodynamic medium at the 

receiver outlet (°F) ..... 

Intercept Factor 

Conduction loss (kW t /m 2 ) . 

Convection loss (kW^/m 2 > . . . . 

Reradiation loss (kWt/m 2 ) ..... 

Combined losses (kW t /m 2 ) 

Efficiency . . . 

Annual Operating Parameters 

Efficiency ..... 

Thermal Transport 

Design Efficiency 

Annual Efficiency 

Power Conversion Efficiency 
Design Parameters 

Heat Engine Efficiency . 

Generator Efficiency 

Power Conversion System Efficiency 

Annual Power Conversion 

System Efficiency . 

Power Conditioning and Transport Efficiency 

1 - (Fractional Parasitic Losses) 

System Performance 
Design Parameters 

Receiver Output Thermal Power (kWf.) 

Electric Power Output (kW e ) . 

Efficiency 

Annual Performance 

Efficiency . 

Notes are presented on the following page. 


Performance 

Design Point 

Annual 

li) 



(b) 


isl 

ID 


<e>._ 


(f) 



la 1 

m 



m 




a) 
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Table 1. (Cont'd): Notes on the Performance Summary Input Data Form 


(a) The power conversion assembly package and its supporting struts block 
insolation that would otherwise fall on the concentrator. This blocking 
reduces the performance of the concentrator. 

(b) Between cleanings, the design reflectivity of the concentrator mirrors 
is degraded by the gradual accumulation of dust on the first surface of 
the mirror panels but is virtually completely restored by proper 
cleaning. 

(c) The annual performance of the usual optimal concentrator field is 
degraded by the shading of some concentrators by others at times 
considerably removed from solar noon. This is particularly true during 
the winter. 

(d) The apparent absorptance for concentrated insolation of the receiver 
apertures considered here is assumed to differ negligibly from unity. 

(e) The maximum bulk temperature achieved by the thermodynamic medium as it 
moves through the complete cycle is noted here. 

(f) Thermal energy is lost from the receiver through the parallel paths of 
( l) reradiation through the aperture, (2) direct convection at the 
aperture, and (3) conduction through the body and/or supports of the 
receiver followed by convection from the surface of the receiver body. 
One minus the sum of these losses (expressed as fractions of energy 
entering the receiver cavity) represents the transfer factor for energy 
entering the receiver aperture. 

(g) Receiver annual efficiency is affected by the insolation histogram 
employed, even in the case of a non-storage system. 

(h) Because the receiver and thermodynamic engine subsystems are closely 
coupled in a point-focusing thermodynamic power module, the module's 
thermal transport efficiency differs negligibly from unity. Annual 
thermal transport efficiency is assumed equal to the design thermal 
transport efficiency. 

(i) Power conversion annual efficiency is affected by the insolation 
histogram employed, even in the case of a non-storage system. 

(j) Annual system efficiency is a function of both receiver and power 
conversion annual efficiencies and thus is affected by the insolation 
histogram employed, even in the case of a non-storage system. 


SECTION 1. 1 


DESCRIPTION OF THE PROBLEM AND ITS SOLUTION 


A. INPUT 

VARIABLES 

(All in Segment No. 1) 

Algebraic 

Variable 

Acronym 

Definition 

X DND 

IDND 

Deaign level of direct normal insolation (often 
assumed to be 0.8 or 1.0 kW/m^). 

^CD 

BCD 

Design efficiency of the concentrator. 

Vc 

HC 

Annual average efficiency of the concentrator 
(degraded by mirror soiling and by the mutual shading 
of grouped concentrators). 

Qrl 

QRL 

Receiver combined thermal loss in kW c /m2 of 
concentrator ares projected onto the concentrator 
aperture plane. 


FSBS 

Intercept factor of the receiver aperture, i.e., the 
fraction of concentrated solar radiation falling on 
the aperture plane of the receiver that passes through 
the aperture. 

a 

ALPH 

Apparent absorptivity of the receiver aperture for the 
concentrated insolation incident on it. 

V PD 

HPD 

Design efficiency of the engine. 

a n 

ANPS 

The degree of the polynomial representation must be no 
greater than 6. 


NOTE: 

PPX-59 Professional Program 208008 Polynomial 
Regression by T. H. Wysmuller has been most useful in 
evaluating the required polynomial coefficients from 
graphs of normalized power conversion efficiency as a 
function of normalized input thermal power* 

Vegt 

HEGT 

Efficiency of the combined electrical generator, and 
transport and power conditioning subsystems* 

^P 

HPAR 

Correction factor by which the gross output of the 
pov;er plant must be multiplied in order to account for 
losses due to plant auxiliary equipment, office and 
control room space conditioning, etc. 
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Algebraic 


Definition 

Variable 

Acronym 

*DN 

IDN 

Actual level of direct normal insolation at which 
performance is to be determined (kW/m^). 

H 

urs 

Annual total hours of direct normal insolation at 
intensity centered at Ipfj CAIqn d 0*05 kW/m^), 

N 

NPOL 

Degree of the polynomial representing the normalized 
power conversion subsystem efficiency. 

0 < N < 6 (Integer). 


B. OUTPUT VARIABLES 


L. In Segment No. 1 


Algebraic 

Variable 

Acronym 

Definition 

^CD 

HOD 

Concentrator design efficiency (an input variable). 

%) 

HRD 

Receiver design efficiency. 

^PD 

HPD 

Power conversion subsystem design efficiency defined 
as thermal input to electrical generator output (an 
input variable) . 

’to? 

HEP 

Efficiency of electricui power conditioning and 
transport including an allowance for module or plant 
parasitic power requirements. 

^SD 

HSD 

System design efficiency. 

2. In 

Segment No. 

2 

Algebraic 

Variable 

Acronym 

Definition 


At each 

value of IDN: 

^ C 

HC 

Concentrator efficiency. 

’JCN 

HCN 

Normalized concentrator efficiency, llC/HCD. 

Qc 

QC 

Concentrator photon power output per unit concentrator 
area (kW/rn^'). 

Qcn 

QCN 

Normalized concentrator photon power output, QC/QCD. 
(See page 2-6 for the definition of QCD. ) 

^R 

HR 

Receiver efficiency. 
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Algebraic 


Variable 

Ac ronym 

Vm 

URN 

q r 

QR 

Qrn 

QRN 

v? 

HP 

Vl'N 

1IPN 

Qp 

QP 

QpN 

QPN 

Vs 

na 

VSH 

HSN 

Qs 

QS 

Qsn 

QSN 


Annual I 

*?RA 

HRA 

^PA 

HPA 

2e l 

2ei 

Vsk 

HSA 


Definition 

Normalized receiver efficiency, IIR/HRD. 

Receiver thermal power output per unit concentrator 
area (kW/m 2 ). 

Normalized receiver thermal power output, QR/QRD. 

Power conversion subsystem efficiency. 

Normalized power conversion subsystem efficiency, 
UP/HPD. 

Power conversion subsystem electrical power output per 
unit concentrator area (kW/m 2 ). 

Normalized power conversion subsystem electrical power 
output, QP/QPD. 

System efficiency. 

Normalized system efficiency, HS/HSD, 

System electrical power output per unit concentrator 
area (kW/m 2 ). 

Normalized system electrical power output, QDS/QSD, 


Annual average receiver efficiency. 

Annual average power conversion subsystem efficiency. 
Annual direct normal insolation [kW/(m 2 y)]« 

Annual average system efficiency. 


C. METHOD OF SOLUTION 

As indicated in both the abstract and the introduction, the program 
presented here has been developed to allow calculation of the design-point and 
quasi-steady-state annual performance of a modular solar thermal electric 
power plant, each module of which comprises (1) a sun-tracking, point- 
focusing concentrator and (2) a power conversion assembly that is mounted on 
the concentrator with the receiver aperture near the focal plane and on the 
optical axis of the concentrator. In turn, the power conversion assembly 
comprises a cavity receiver, a thermodynamic engine complete with auxiliaries, 
and an electrical generator. The electrical control and conditioning 
equipment directly associated with each modular power conversion assembly may 
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be concentrator-mounted or ground-mounted, or ports of it may be 
concentrator-mounted and the remainder, ground-mounted. Typically, the power 
plant in its entirety will require additional power control and conditioning 
equipment to enable connection to a power grid and to supply auxiliary power 
to the modules, at least during start-up. 

It is again emphasized that the program developed here is valid only for 
a power plant that has negligible energy storage capacity for either thermal 
or electric energy. 

Another more subtle but ordinarily less important limitation of the 
program is that the plant power output may exceed its design value if the 
direct normal insolation exceeds the design value selected for this 
parameter. MorG sophisticated programs ordinarily impose the requirement that 
the design output of the power plant never be exceeded, wasting the excess 
power if necessary. However, because of TI-59 limitations in available 
program steps and memory cells, the author has been unable to incorporate this 
feature into the program presented here. If the direct normal insolation 
received at a given site seldom exceeds the selected design value or only 
exceeds it by small amounts, this limitation on the program is negligible. On 
the other hand, if a design direct normal insolation of 800 W/m 2 were chosen 
for a plant whose power conversion subsystem could accept the energy collected 
by the concentrator at a direct normal insolation of 1000 W/m 2 and if this 
plant were located in the desert Southwest, the difference in annual 
performance as calculated by this program and by one limiting output to the 
design maximum would be appreciable. 

Figure 1 has been included so as to provide a graphic representation of 
the energy path through the power plant. Most of the direct normal insolation 
collected, reflected, and concentrated by the concentrator is directed to the 
receiver os photon energy; most of the photon energy absorbed by the receiver 
is transferred as thermal energy to the thermodynamic engine; the 
thermodynamic engine converts some of the thermal energy to mechanical energy 
and rejects or loses the remainder; the generator converts most of the 
mechanical energy to electrical energy with only slight energy dissipation; 
the module's electrical power control and conditioning transfers the energy to 
the plant control and conditioning system, again with relatively small losses; 
finally, the plant control and conditioning system transfers most of this 
energy to the electrical grid, while most of the remainder is absorbed by 
plant auxiliaries and the rest is dissipated. 

Concentrator performance is determined in the following manner. At 
design conditions, 


a cQcd " A C I DND/ > C F B in kw 


represents the photon flux transmitted to the receiver aperture plane by the 
concentrator at design conditions, where Aq represents the concentrator's 
reflecting surface area projected on the plane of the concentrator aperture, 
Pq represents mirror reflectance, and F3 represents the fraction of Aq 
that is blocked by the power conversion assembly package and its supporting 
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struts. Ip ND is defined in Section A of this chapter. Under general 
conditions , 

A cQc “ a c i dn/ 5 c f b • f dg ’ F s * 


where Fn G represents the average fractional degradation of mirror 
reflectivity that results from dust accumulation on the first surfaces of the 
mirrors between cleanings, and Fg represents the average degradation of 
concentrator field performance for the entire power plant due to the shading 
of some concentrators by others. Even in an optimally distributed field of 
concentrators, shading is appreciable near sunrise and sunset and, at 
latitudes appreciably above or below the equator, is more important in the 
winter than in the summer. Obvious definitions can be written for the design 
and general efficiencies of the concentrator as 


^CD e ^C f B 


and 


1c ~ / j C F B F DG F S > 


respectively. Then Qqd and can be expressed by 


Qcd ° 7 ?cd i dnd 


and 


Qc " ^C^DN * 


Note that Q CD , Qq, and all other Q's in the following discussion represent 
energy flux per unit area of concentrator reflective surface projected upon 
the aperture plane of the concentrator. Defining the Q's in this way allows 
any one of them along the path of energy transport and conversion to be 
compared directly with the incident direct normal radiation, a 

normalized photon flux from the concentrator, valid for annual average 
conditions of shading and degradation, can then be defined for an arbitrary 

X DN 3 3 


Qcn ° Qc/Qcd • 
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The definition of receiver efficiency employed by the program is more 
complicated than that for the concentrator. Receiver thermal flux delivered 
to the thermodynamic engine is defined in the following wayi 


Qrd = g<£Qcd “ Qrl > 


where all the parameters and variables on the right side of this equation have 
already been defined. Defining receiver design efficiency as 

„ Q RD 


by substitution, 


« ~ nA, RL 


At other than design conditions, the receiver thermal flux delivered to 
the engine is defined as 


Qr ~ a< ^Qc “ Qrl > 


and receiver design efficiency can be defined as 


^R " Q r 


= a</> - 


q rl 

Qc 


Then from earlier definitions, 


^RD ~ ^R " - Q 


X RL 

'CD 


1 “ 


q c/ q cdJ 


and 


Qrl 

^R^RD + Q^ 


1 - 


Qf! 


CNy 
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Finally, a normalized receiver efficiency can be defined as 


’JRN “ '’W^RD * 

As can be seen from the formulae for r/%, at insolation levels below a 
certain threshold, receiver thermal efficiency can become negative. The 
program presented here sets 17^ 0 0 whenever 17^ < 0 is calculated; 
obviously, all subsequent Q's and tj's except Tjp are set equal to zero for 
insolation levels below this threshold. 

The evaluation of the energy flux through the remainder of the power 
conversion subassembly is treated in a manner different from that employed in 
the case of the concentrator or receiver. The full-load efficiency of the 
engine is considered known from extensive analysis and/or experiment by the 
manufacturer. The same knowledge is usually available for the generator and 
the complete power conditioning and transport system. The fractional 
parasitic losses for the plant are also known from experiment or have been 
estimated. 

In the program presented here, parasitic losses are assumed to be a 
small constant fraction of the gross plant-produced, conditioned, and 
transported electrical output; thus, a constant component efficiency 
reflecting plant parasitic losses can be defined as 


’’/PAR “ 1 " f PAR • 

where fp AR represents the normalized fractional parasitic losses just 
described. The efficiencies of the generator and of the electrical power 
conditioning and transport subsystem can be entered during the first segment 
of program execution as normalized fractional constants. Typically, real 
generator and electrical conditioning and transport subsystems exhibit 
efficiencies that reflect a relatively small, constant, fractional no-load 
loss and an additive fractional loss that varies with the fraction load. This 
same behavior is also exhibited by the typical thermodynamic engine. 

Therefore, if the user of this program wishes a more accurate representation 
of the combined efficiency of the engine, generator, and power conditioning 
and transport subsystems and if combined efficiency as a function of load is 
known or can be estimated for this combination, a polynomial representation 
for this efficiency, Tjp ♦ ^ 7 eGT> can entered during the loading 
phase (Segment No. 1 ) of the program in place of 7 )^ - and ^EGT ~ ^ 

can be set. Another alternative often employed by the author is to define the 
normalized efficiency of the engine-generator combination by a least-squares- 
best-fit polynomial function of normalized thermal flux from the receiver to 
the engine (converted from the manuf actu7:er 1 s data presenting efficiency as a 
function of output for this combination) and to assume respective constant 
values for 7 and fp^R> the efficiency of the power conditioning and 
transport subsystem, and the fractional parasitic losses, respectively. 
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In the sample problem and corresponding solution here, T7EGT and ''PpAR 
have been defined as constants, and only the efficiency of the engine has been 
defined as a function of the heat flow into the engine. 

The program presented here allows the representation of normalized 
engine efficiency as a polynomial function of normalized thermal flux to the 
engine from the receiver. This representation is limited to a polynomial of 
no more than sixth degree. 

Once the normalized thermal flux from the receiver, Qrn, has been 
determined, the normalized efficiency of the engine, ^p^, can be 
calculated. Then the normalized mechanical output power of the engine, its 
efficiency, and its dimensional mechanical output power can be calculated: 


Qpn “ ^PN * Qrn * 

’JP " ^PN * ^PD » 
Qp = t? p • Qr . 


As was the case with the receiver, at input thermal fluxes to the engine 
below some threshold, 17 p will be negative. When this occurs, the program 
sets 7 jp^ = 0, and all subsequent Q's and 77 ' s in the energy train are also 
set equal to zero. 

If the product of 17 egt efficiency of the generator and electrical 

power conditioning and transport subsystem) and (the normalised 

fractional correction for plant parasitic electrical power dissipation) is 
defined as ^gp, the electrical power output of the plant can be calculated 
for design and arbitrary loads as 


Qsd “ ^EP ’ Qpd 

and 

Qs “ ^EP * Qp • 


For design and arbitrary loads, plant efficiency can be obtained from the 
simple formulae, 


Isd = Qsd^dnd 


and 


Vs ~ Qs^dn • 
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A formula for Che normalized Bystem efficiency of Che plant is 


^SN 13 Vs/VSD • 


The assumption of quasi-steady plant operation allows all transients due to 
abrupt changes in insolation to be ignored; thus, it becomes unnecessary to 
employ a dynamic temporal record of direct normal insolation at a particular 
site in order accurately to simulate annual plant performance. A histogram of 
annual direct normal insolation at the same site serves just os well under 
this assumption, and, instead of making thousands of calculations of plant 
performance to obtain annual plant performance, only a few tens of 
calculations of plant performance are required. Histograms for a few typical 
sites during what have been judged to be representative years are presented in 
Figure 2. 

Annual plant performance for a plant is calculated by the program in the 
following way. At the average IpN corresponding to each incremental Alpfj , 
energy fluxes are calculated and multiplied by the number of hours during the 
year in which that particular was measured* 

E 1 “ IDN * » 
s c “ Qc * H 

E R = Q R . H 
Ep » Qp • H 

E s = Qs * H 


These individual energies are then added to the sum of corresponding energies 
determined for all the previously employed (smaller) values of average Ipjj. 
When all bands of I D ^ have been considered, summations for the year are then 
available : 


SEj “ 2 (Iqn • H) 
2E c = Z(Q C • H) 
Ze r = 2(Q r . h) 
SE p = S(Q p • H) 
ZE S = Z(Q S . H) 
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Figure 2. Histograms of the Annual Hourly Average Direct Insolation 
for Four Selected Sites 
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Finally, from these annual summations, the annual average efficiencies of the 
receiver, the engine, and the system can be calculated: 


^RA “ 215 r /Ze c 

^PA B ^Ep/ZE R 
Vsh a ^Eg/SEj 


The annual average efficiency of the concentrator was established near the 
beginning of this section. It is therefore unnecessary to determine 


^CA “ ^Eq/SEj , 


although this equation can be employed as a check. 


D. FLOW CHART 

The flow chart (Appendix A) of the programmed solution follows standard 
flow charting conventions and is virtually self-explanatory. Segment No. 1 is 
stored on the first TI-59 program card and Segment No. 2, on the second. The 
names of the various subprograms correspond to the various TI-59 keys employed 
as user-defined subprogram labels in the course of program development. 
Although an attempt was made to choose each label in a way that would provide 
a hint as to the function of the corresponding subroutine, the severely 
limited number of names available probably prevented much success in this 
respect. The reader must simply remember the names of the various subprograms 
as the various loops of the program are traced and retraced. 


E. ANNOTATED PROGRAM LISTING 

While the flow chart presented in Appendix A is easily adapted to any 
programmable calculator or microcomputer whose capabilities equal or exceed 
those of the TI-59, the program listing presented in Appendix B can be 
directly keyed in and executed only on a TI-59. Nevertheless, since a TI-59 
and attached printer are still commonly available and widely used in science 
and engineering, the program listing has been included in this report. To 
those who are familiar with the TI-59 programmable calculator and the 
ancillary, no-longer-extant PPX-59 user's group, the form in which the listing 
is presented will be completely familiar. The comments accompanying the 
listing closely parallel the remarks presented in the flow chart. 


SECTION III 


SAMPLE PROBLEM AND SOLUTION 


The design and annual performances of a paraboloidal concentrator 
(dish)/Brayton solar thermal power plant (PDB/STPP) operating without thermal 
or electrical storage capability at Barstow, California, is considered here. 

A histogram of the direct normal insolation received at Burstow in 1976 is 
presented in tabular form as Table 2. 

The design value for direct normal insolation is assumed to be 
1.0 kW/m^. Design and annual values of concentrator efficiency are assumed 
to be 0-.94 and 0.88, respectively. 

A receiver outlet temperature of 870°C (1600°F) has been specified. The 
corresponding combined thermal loss, intercept factor, and apparent aperture 
absorptance for the receiver have been estimated to be 0.07445 kW/m^ (based 
on concentrator aperture), 0.99, and 0.92, respectively. The design 
efficiency of the Brayton engine has been estimated as 0.271, and the 
normalized part-load characteristics of the engine are presented in Table 3. 


Table 2. Tabular Histogram of Direct Normal Insolation Received 
at Barstow, California, in Calendar Year 1976 


Insolation Hours of Insolation 

Increment in this Increment 


0.00 

- 

0.05 

251 

0.05 

- 

0.10 

174 

0.10 

- 

0.15 

131 

0.15 

- 

0.20 

109 

0.20 


0.25 

111 

0.25 

- 

0.30 

122 

0.30 

- 

0.35 

94 

0.35 

- 

0.40 

92 

0.40 

- 

0.45 

96 

0.45 


0.50 

111 

0.50 

- 

0,55 

141 

0.55 

- 

0.60 

129 

0.60 

- 

0.65 

143 

0.65 

- 

0.70 

192 

0.70 

- 

0.75 

228 

0.75 

- 

0.80 

288 

0.80 

- 

0.85 

345 

0.85 

- 

0.90 

489 

0.90 

- 

0.95 

595 

0.95 

- 

1.00 

508 

1.00 

- 

1.05 

60 


3-1 



Tabic 3. Estimated Garrett Turbine Engine Company SAGT-1A 
Drayton Engine Normalized Pert-Load Efficiency 
As a Function of Normalized Input Thermal Power 


Normalized 
Input Thermal Power 

Normalized 

Part-Load Efficiency 

0.295 

0.000 

0.375 

0.515 

0.500 

0.827 

0.625 

0.940 

0.750 

0.989 

0.875 

1.000 

1.000 

1.000 


The coefficients of the sixth-degree polynomial approximation 
corresponding to this normalized part-load operating characteristic are as 
follows : 

a 0 >=> -11.484 
a! “ +93.72,1 

a 2 = -303.71 
a 3 n +532.70 
a^ * -526.26 
a^ = +257.56 
a 6 « -59.527 


The combined efficiency of mechanical-to-electrical power conversion by 
the generator, the power conditioning equipment, and the electrical transport 
network has been estimated to be 0.95, while the correction factor to the 
gross output of the power plant required to account for losses due to plant 
auxiliary equipment, etc., has been estimated to be 0.98. 

As indicated by the sample problem (Appendix C), the input segment of 
the program requests the entry of the input data, item by item. Entering the 
degree of the polynomial (NPOL) and starting the program results in 
(1) calculation and printing of the design performance parameters 
corresponding to the PDB/STPP configuration specified by the input parameters 
and (2) a printed directive to load the second segment of the program into the 
calculator. 
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After the second segment is loaded and started, the annual performance 
of tho specified PDB/STPP is calculated. When the calculator stops, the 
,lution is complete. The following results have then been presented: 

H) performance data for the concentrator, receiver, engine, and complete 
system for each median value of direct normal insolation employed in the 
histogram — only the results obtained at In^ ° 0.025 are labeled — and 
(2) the total direct normal insolation received during the year par square 
meter of concentrator aperture area and the annual average efficiencies of the 
receiver, thermodynamic engine, and the system in its entirety. 
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SECTION IV 


USER INSTRUCTIONS 


While the instructions presented in Appendix D arc intended primarily 
for instruction of the TI-59 user and while the contents of the data registers 
are uniquely applicable to the TI-59, the general sequence of instructions Ib 
certainly applicable to any programmable calculator or microcomputer that 
might be employed to solve the problem* 
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SECTION V 


CONCLUDING REMARKS 

Tho results obtained for a wide variety of paraboloidal-dish solar 
thermal power plnnto operating at vnriouo oiteo have been compared with 
hour-by-hour calculations made for the same configurations at tho same sites. 
In every case for which the site was located in tho sunbelt and for which 
*DND ** was specified, tho results obtained by the method described 

herein differed negligibly from the hour-by-hour resultB obtained with a 
mainframe computer. 

For the reader who has access to a TI-59 programmable calculator and 
attached TI PC-100C printer, loading the program and duplicating the sample 
problem should be a relatively simple task. 

For the reader who wishes to understand the problem and its solution, 
the task i3 more difficult. However, the author recently presented the flow 
chart and the TI-59 program listing that are included here to members of a 
senior class in mechanical engineering. Within a few days these students had 
rewritten the program in BASIC and within tan days had duplicated the sample 
problem results with a Commodore 64 computer. Since that time, these same 
students have successfully programmed and solved the sample problem (1) on an 
HP-41C programmable calculator and (2), using FORTRAN, on the university's 
mainframe computer. Thus, the author is convinced that, by studying this 
report, un interested reader who is familiar with a sufficiently capable 
programmable calculator or computer will be able to create a program and solve 
the sample problem presented herein and any similar problems of interest. 
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APPENDIX A 


FLOW CHART OF THE PROGRAMMED SOLUTION 


A-l 









DR 51 > DSPR 

Note: DSPR designates the 
display register 

Print DSPR 

Note: DR51 = HRN if Flag 0 
has not been set. 

DR 51 * HPN if Flag 0 
has been set. 
















































APPENDIX B 


TI-59 PROGRAMMABLE CALCULATOR PROGRAM 
SEGMENT NOS. 1 AND 2 


B-l 


Segment No. } 


JRIGINAL P^v. ' 
OF POOR QUAUrir! 


IOC CODE KEY 


COMMENTS 


LOC 

CODE 

KEY 

055 

90 

LST 

056 

01 

1 

057 

03 

3 


COMMENTS 


LOC CODE KEY 


COMMENTS 


FSBS +DR39 
A 


110 36 36 


HCD 


000 

03 

o 

o 


001 

04 

4 


002 

42 

STO 


003 

00 

00 

34+ DROO 

004 

71 

SBR 


005 

99 

PRT 


006 

02 

2 

I 

007 

04 

4 

003 

01 

1 

n 

009 

06 

6 


010 

03 

3 

N 

Oil 

01 

1 

012 

01 

1 

n 

013 

06 

6 

u 

014 

71 

SBR 

Print 

015 

90 


ENTER IDND 

016 

02 


IDND + DR35 

017 

03 


H 

018 

01 


r 

019 

05 


b 

020 

01 


D 

021 

06 


022 

00 



023 

00 


Print 

024 

b* 

71 

90 


ENTER HCD 
HCD + DR36 

026 

02 

V 

u 

H 

027 

03 

3 

023 

01 

i 

r 

029 

05 

5 

V 

030 

00 

0 


031 

00 

0 


032 

00 

Ci 


0 3 8 

00 

0 

Print 

034 

71 

SBR 

ENTER HC 

035 

90 

LST 

HC + DR37 

036 

03 

3 

Q 

037 

04 

4 

033 

03 

0 

sJ 

R 

039 

05 

5 

040 

02 

S' 

iL 

L 

041 

07 

7 


042 

00 

0 


043 

00 

0 

Print 

044 

71 

SBR 

ENTER QRL 

0 4 5 

90 

LST 

QRL + DR38 

046 

02 


F 

047 

01 

T 

048 

03 

3 


049 

06 

6 

b 

050 

01 

1 

B 

051 

04 

4 

052 

03 

3 

5 

053 

06 

6 

Print 

054 

71 


ENTER FSBS 


058 

02 

2 

I 

059 

07 

7 

L 

060 

03 

3 

p 

061 

03 

3 

062 

02 

2 

H 

063 

03 

3 

Print 

064 

71 

SBR 

ENTER ALPH 

065 

90 

LST 

ALPH+DR40 

066 

02 

2 

H 

067 

03 

3 

068 

03 

3 


069 

03 

'2 

P 

070 

01 

1 

n 

071 

06 

6 

u 

072 

00 

0 


073 

00 

0 

Print 

074 

71 

SBR 

ENTER HPD 

075 

90 

LST 

HPD + DR41 

076 

02 

2 

H 

077 

03 


078 

01 

1 

r 

079 

07 

7 

t 

080 

02 


G 

081 

02 

V 

u 

082 

03 

o 

T 

0 8 3 
084 

07 

71 

7 

SBR 

i 

Print 

ENTER HPAR 

085 

90 

LST 

HPAR +DR42 

U 8 6 

02 

2 

H 

087 

03 


088 

03 

3 


089 

03 

3 

P 

030 

01 

i 

A 

091 

03 

3 

A 

092 

03 

3 

R 

093 

05 

5 

094 

095 

69 

02 

DP 

02 

Print 

EN1ER HPAR 

096 

69 

□ P 

097 

05 

05 


098 

91 

R/S 

Enter HPAR 

099 

99 

PRT 

Print HPAR 

100 

49 

PRD 

HPAR * 

101 

42 

42 

HEGT+DR42 

102 

43 

RCL 

103 

40 

40 

ALPH * 

104 

49 

PRD 

FSBS+DR39 

105 

39 

39 


106 

43 

RCL 

IDND 

107 

35 

35 


108 

65 

X 

* 

109 

43 

RCL 



111 

95 

— 

= QCD 

1 12 

42 

STO 

113 

40 

40 

3CD+DR40 

114 

98 

RDV 


115 

04 

4 


116 

42 

STO 


1 17 

00 

00 

4+ DROO 

118 

05 

5 


119 

42 

STO 


120 

33 

33 

5+DR33 

121 

01 

1 


1 *>*> 

1 W fart 

01 

1 


1 23 

42 

STD 


124 

34 

34 

11 + DR34 

125 

01 

i 

A 

126 

03 

3 

A 

127 

03 

3 

M 

1 28 

01 

1 

IN 

129 

03 

3 


130 

03 

3 

p 

131 

03 

3 


132 

06 

6 

S 

133 

69 

OP 

Print 

134 

02 

02 

ENTER ANPS 

135 

69 

OP 


136 

05 

05 


137 

76 

LBL 


1 T 1 C r 
JL v 1 0 

78 

1 + 


139 

01 

1 

DROO + 1 

140 

44 

SUM 

+DROO 

141 

00 

00 


142 

43 

RCL 

DROO 

143 

on 

00 

144 

75 

- 

“ 

145 

43 

RCL 


146 

33 

33 

DR33 

147 

95 

= 

=n (of a n ) 

148 

39 

PRT 

Print n 

149 

91 

R/S 

Enter a n 

150 

151 

152 

1 53 

72 

00 

93 

43 

ST* 

00 

PRT 

RCL 

Store a n 
Print a n 
Continue 
entering 

154 

155 

156 

43 

84 
XI T 
RCL 

and storinc 
a n through’ 
n - 6. 

157 

00 

00 

158 


I MV 

Note: If 

159 

77 

GE 

NPOL=m < 6 

160 

78 

2 + 

store a n =0 

161 

i"l i“i 

7 0 

ADV 

for n > m. 

162 

01 

1 


1 6 3 

02 

2 


1 6*4 

42 

STO 
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PRECEDING PAGE ELANK NOT FILMED 




















LOC CODE 


COMMENTS 


12+DR33 


32 + DR34 
H 



Print 
ENTER 
HOURS 
AT 
I DM 


|| LOC |CODE| 

KEY 

COMMENTS 

2 2 Cl 

04 

.04 


221 

01 

1 


v y y 

u. t— 

44 

SUM 

DROO + 1 

223 

00 

00 

•*DR00 

224 

43 

RCL 


225 

00 

00 


226 

75 

- 


227 

43 

RCL 


228 

33 

'*i *“< 
."i . 


228 

95 


Calculate 

230 

65 

X 

IDN 

231 

93 



<-< -~i 

00 

0 


233 

05 

5 


234 

QC 
O V 

+ 


235 

93 



236 

00 

0 


237 

02 



£38 

05 

5 


239 

95 


T 

240 

69 

□ P 


241 

06 

06 

Print IDN 

242 

69 

□ P 


243 

00 

00 


244 

06 

6 


245 

04 

4 


246 

02 

2 

if 

247 

03 

3 

H 

248 

03 

o 

D 

249 

05 

5 

K 

250 

03 

O 

c 

251 

06 

6 

j 

252 

69 

□p 


253 

04 

04 


254 

25 

CLR 


255 

91 

R/S 

Enter HRS 

256 

69 

□ P 

Print HRS 

257 

06 

06 


258 

72 

ST* 

Store HRS 

259 

on 

00 


260 

43 

RCL 

Continue 

261 

34 

34 

storing 

262 

32 

XT T 

HRS 

263 

43 

RCL 

corres- 

264 

00 

00 

ponding to 

265 

22 

I MV 

IDN until 

266 

77 

GE 

the 21 

267 

60 

DEG 

bins have 

268 

98 

fmv 

been filled 

269 

71 

SBR 


270 

99 

PRT 


271 

03 

3 


272 

01 

1 

N 

273 

03 

3 

p 

274 

03 

3 

r 


LOC CODE] KEY 


275 0 


COMMENTS 


HRD+DR43 


oimimi Pmz & 

°P POOR QUALITY 























OF PGC;< C U. 
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Seament N 


LOC |i 

CODE) 

KEY 

ODD 

01 

1 

001 

01 

1 

002 

42 

STD 

003 

00 

00 

004 

76 

L SL- 

005 

69 

OP 

006 

22 

INV 

007 

86 

STF 

003 

00 

00 

009 

98 

RDV 

010 

69 

DP 

Oil 

00 

00 

012 

06 

6 

013 

04 

4 

014 

02 

'■t 

i»< 

015 

04 

4 

016 

01 

1 

017 

06 

6 

018 

03 

3 

019 

01 

1 

020 

69 

DP 

021 

04 

04 

022 

01 

1 

023 

44 

SUM 

024 

00 

00 

025 

QO 

r o 

» 

026 

00 

o 

027 

02 

2 

028 

05 

5 

029 

85 

+ 

030 

Cl O 

a 

031 

00 

0 

032 

05 

nj 

033 

65 

X 

034 

53 

< 

035 

43 

RCL 

036 

00 

00 

037 

75 

- 

038 

01 

1 

039 

02 

2 

040 

54 

) 

041 

95 

= 

042 

69 

OP 

043 

06 

06 

044 

98 

RDV 

045 

42 

STD 

046 

46 

46 

047 

43 

RCL 

048 

37 

37 

049 

99 

F'RT 

050 

55 

■7 

051 

43 

RCL 

052 

o c 

Q 

O £ 
O D 

053 

95 

— 

054 

99 

F’RT 


COMMENTS 


U+DRQQ 


Clear 

Flag 

0 


DROO + 1 
-►DROO 


Calculate 

IDN 


Print 

X.XXX=IDN 

IDN+DR46 


Print HC 
HC 


=HCN 

Print HCN 


LOC CODE 

KEY 

COMMENTS 

LOC CODE 

KEY 

055 

43 

RCL 


1 10 

43 

RCL 

056 

46 

46 

IDN 

111 

44 

44 

057 

65 

X 

* 

112 

95 

— 

058 

43 

RCL 


113 

99 

PRT 

059 

37 

37 

HC 

114 

98 

RDV 

060 

95 

zz 

=qc 

115 

76 

LBL 

061 

99 

PRT 

Print QC 

116 

44 

SUM 

062 

42 

STD 

QC-t- DR47 

117 

42 

STD 

063 

47 

47 

QC 

118 

03 

03 

064 

55 

•7 

• 

119 

53 

( 

065 

43 

RCL 

QCD 

120 

43 

RCL 

066 

40 

40 

8 QCN 

121 

04 

04 

067 

95 

— 

Print QCM 

122 

42 

STD 

068 

99 

PRT 

QCN -*-DR4Q 

123 

02 

02 

069 

42 

STD 

124 

85 

+ 

070 

48 

48 


125 

05 

5 

071 

98 

RDV 


126 

54 

> 

072 

35 

l/H 

1/QCN 

127 

42 

STD 

073 

75 

- 


128 

01 

01 

074 

01 

1 

1 

129 

01 

1 

075 

95 

S2 

c 

130 

44 

SUM 

076 

65 


(1/QCN-l) 

131 

02 

02 

077 

43 

RCL 

* 

132 

73 

RC* 

078 

w* O 

o o 

• Z' o 

QRL 

133 

01 

01 

079 

55 

•f 


134 

^ 6 

LBL 

080 

43 

RCL 

QCD 

135 

60 

DEG 

081 

40 

40 


1 36 

y y 

INV 

082 

75 

- 


137 

97 

DSZ 

083 

43 

RCL 

HRD 

138 

02 

02 

084 

43 

43 

— UD 

139 

70 

F'RD 

085 

95 

= 

- - HK 

140 

69 

DP 

086 

94 

+/ - 

HR 

141 

31 

31 

087 

99 

PRT 

Print HR 

142 

53 

< 

088 

42 

STD 

HR + DR49 

143 

24 

CE 

089 

49 

49 

HR 

144 

65 

v 

r*. 

090 

55 

•f 

• 

145 

43 

RCL 

091 

43 

RCL 

HRD 

146 

03 

03 

092 

43 

43 

=HRN 

147 

85 

+ 

093 

95 


148 

73 

RC* 

094 

42 

STD 

HRN-*-DR51 

149 

01 

01 

095 

096 

51 

61 

51 

GTD 

Go to 

150 

151 

54 

61 

> 

GTD 

097 

48 

EXC 

Subroutine 

152 

60 

DEG 

098 

76 

LBL 

EXC 

153 

76 

LBL 

099 

49 

F'RD 


154 

70 

RRD 

100 

43 

RCL 


155 

42 

STD 

101 

47 

47 

QC 

156 

51 

51 

102 

65 

X 

* 

157 

65 

X 

103 

43 

RCL 

HR 

158 

43 

RCL 

104 

49 

49 


159 

41 

41 

105 

95 

— 

-QR 

160 

95 

zz 

106 

99 

PRT 

Print QR 

161 

99 

PRT 

107 

42 

STD 

QR-» DR50 

162 

42 

STD 

1 08 

50 

50 

QR 

163 

52 

52 

109 

55 

■7 

• 

164 

86 

STF 


COMMENTS 

QRD 

=QRN 

Print QRN 


QRN + DR03 


NPOL 

NPOL DR02 


NPOL + 5 
*DR01 
DR02 + 1 
+ DR02 

DR(DROl) 

8 ag + DSPR 

If DR02 n 0 
go to 

Subroutine 
RAD; other- 
wise, con- 
tinue. 

DR01 - 1 
-*• DROl 

DSPR 

* 

QRN 

+ 

DR (DROl) 


Go to 

Subroutine 

DEG 

HPN+DR51 i 
HPN 
* 

HPD 

=HP 

Print HP 
HP + DR52 
Set 


















ORIGINAL PAGE *9 
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LOC (CODE 

KEY 

COMMENTS 

|| LOC CODE 

KEY 

COMMENTS 

LOC 

taaii 

KEY 

["comments 1 

165 

00 

00 

Flog 

220 

43 

RCL 


275 

57 

57 


166 

61 

Gin 

0 

221 

03 

03 

Print QPN 

276 

55 

- 

t 

16? 

168 

48 

76 

EMC 

LBL 

Go to 

Subroutine 

222 
km ^ •»' 

99 

43 

PRT 

RCL 

"9 

Ci * 

278 

43 

56 

RCL 

56 

EQC*HRS 

1 69 

59 

I NT 

EXC 

224 

46 

46 


279 

95 

SI 

n HRA 

i?o 

43 

RCL 

QRN 

X 

225 

71 

SBR 


280 

69 

DP 

Print 

171 

03 

03 

226 

39 

ir 

EIDN*HRS 

281 

06 

06 

ZZXX°HRA 

172 

65 


n 

© © 7 

Vm N 1 

44 

SUM 

282 

98 

flDV 


173 

43 

RCL 

HPN 

„nr»M 

© ©•© 
CC* 

55 

55 

B EEI 

km ’J 

69 

DP 


174 

51 

51 

n qPN 

©©q 

Lm km m* 

43 

RCL 

•+-DR55 

284 

00 

00 


175 

95 

= 

now mk, non^j 

230 

47 

47 


285 

06 

6 


176 

42 

STD 

L|rN *r UKU J 

231 

71 

SBR 


286 

04 

4 

s 

177 

03 

03 


232 

89 

IT 


287 

02 

y 

u 

178 

43 

RCL 

HPN 

233 

44 

SUM 

EQC*HRS 

288 

03 

3 

H 

173 

51 

51 

Y, 

234 

56 

56 

+DR56 

289 

03 

3 

n 

180 

65 

it 

©©■=: 
k# ’J 

43 

RCL 


290 

03 

3 

P 

181 

43 

RCL 

HPD 

©3 A 

km 'm> © 

50 

50 


291 

01 

1 

A 

182 

41 

41 

© © 7 

km W ( 1 

71 

SBR 


292 

03 

3 

f\ 

183 

95 

= 

c HP 

233 

89 

if 

EQR*HRS 

293 

69 

DP 


184 

42 

STD 

HP + DR52 

239 

44 

SUM 

294 

04 

04 


185 

52 

52 

HP 

aw 

mm 


+DR57 

295 

43 

PCL 

SQP*HRS 

186 

65 


X 

% 

19 



296 

58 

tr © 
-J © 


1 0 7 

43 

RCL 

QR 

ecrJ 

El 



297 

55 

■r 

*r 

188 

50 

50 

= nn 

243 

71 

SBR 


298 

43 

RCL 

SQR*HRS 

189 

95 

99 

s 


244 

© Cl 
© 7 

tr 

EQP*HRS 

299 

57 

57 

190 

PRT 

Print QP 

245 

44 

SUM 

300 

95 


=HPA 

191 

42 

STD 

tr © 

246 

58 

58 

-►DR58 

301 

69 

DP 

Print 

192 

52 

QP + DR52 

247 

43 

RCL 


302 

06 

06 

ZZXY=HPA 

193 

43 

RCL 

QPN ** 

248 

53 

53 


303 

98 

flDV 


194 

03 

03 

Print QPN 

249 

71 

SBR 


304 

69 

DP 


195 

99 

PRT 

250 

89 

If 


305 

00 

00 


196 

43 

RCL 

QP 

251 

44 

SUM 

EQS*HRS 

306 

06 

6 

B 

197 

52 

52 

252 

59 

era 
■j ” 

307 

04 

4 


198 

65 

X 

X 

253 

03 

3 

“♦■DR59 

308 

07 

7 


199 

43 

RCL 

HEP 

254 

02 

© 

Cm 

If 

309 

07 

7 


200 

42 

42 

255 

32 

MIT 

310 

01 

1 

p 

201 

95 


=QS 

256 

43 

RCL 

DROO < 32 

311 

07 

7 

L* 

202 

203 

98 

42 

flDV 

STD 

QS -»DR53 

nc 

257 

258 

00 

22 

00 

INV 

go to 

Subroutine 

312 
2 1 2 

02 

04 

2 

4 

i 

204 

53 

53 

qs 

259 

77 

GE 

OP; 

314 

69 

DP 


-;05 

55 

*r 

« 

t 

260 

69 

DP 

otherwise » 

315 

04 

04 


206 

43 

RCL 

IDN 

261 

98 

flDV 

contl nue. 

316 

43 

RCL 


20 ? 

46 

46 

e HS 

262 

69 

DP 


317 

e tr 

55 

55 

Print 

208 

95 

52 

263 

00 

00 


318 

69 

DP 

ZZXZ = EEI 

209 

99 

PRT 

Print HS 

264 

06 

6 


319 

06 

06 


210 

42 

STD 

HS+DR54 

265 

04 

4 


320 

98 

flDV 


211 

54 

54 

HS 

266 

02 

2 

H 

321 

69 

DP 


© 1 © 

55 


• 

267 

03 

3 

n 

322 

00 

00 


£ 1 3 

43 

RCL 


268 

03 

.-i 

O 

R 

323 

06 

6 


214 

45 

45 

HSD 

2 69 

05 

Cj 


324 

04 

4 


•“i * cr 

£15 

95 


=HSN 

270 

01 

i 

A 

3 2 5 

02 

o 

u 

u 

216 

99 

PRT 

Print HSN 

271 

03 

3 

n 

326 

03 

3 

H 

217 

43 

RCL 


272 

69 

DP 


327 

03 

3 

s 

218 

53 

53 

nt* i 

273 

04 

04 


328 

06 

6 


219 

99 

PRT 

Print qs 

274 

43 

RCL 


329 

01 

1 
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LOC (CODE 


COMMENTS 


LOC CODE 


COMMENTS 


LOC ICODE 


COMMENTS 


230 

03 

3 

231 

69 

OP 

Vtffa. 

04 

04 

<*'* 11 * 

w* V 

43 

RCL 

234 

59 

59 

335 

55 

t'm 

336 

43 

RCL 




mm 



339 

69 

DP 

240 

06 

06 

241 

91 

R.'S 

342 

76 

LBL 

343 

48 

ehc 

244 

00 

0 

345 

32 

JUT 

846 

43 

RCL 

347 

51 

51 

348 

77 

GE 

349 

87 

IFF 

350 

©-? 
3 i 

IFF 

<251 

00 

00 

•*i l? 

, •! £ 

39 

CDS 

353 

69 

□ P 

354 

00 

00 

355 

02 

2 

356 

33 

3 

5 7 

03 

3 

sT* J O 

05 

5 

859 

03 

3 

360 

01 

1 

361 

06 

6 

362 

04 

4 

363 

00 

0 

364 

01 

1 

365 

69 

DP 

366 

01 

01 

367 

03 

3 

368 

06 

6 

369 

01 

1 

370 

07 

7 

371 

03 

3 

372 

07 

7 

373 

00 

0 

374 

00 

0 

375 

69 

DP 

376 

02 

02 

377 

69 

DP 

378 

05 

05 

379 

00 

0 

380 

42 

STD 

*7 i © 4 

» O X 

49 

49 

© C r o 
O 0 & 

61 

GTD 

383 

49 

PRD 

384 

76 



A 


EQS*HRS 

« 

EEI 

°HSA 

Print 

HZVX°HSA 

End of 
Segment 
No. 2 
If 

DR5U0 
go to 

Subroutine 
IFF; other- 
wise, con- 
tinue. If 
Flag 0 is 
set, go to 
Subroutine 
COS lather- 
'd se, continue 
H 

R - 
N 


0 


S 

E 

T 


Print 

HRN= 0 SET 


0+DR49 
Go to 

Subroutine 

PRD 


385 

386 

387 


1:90 

391 

392 

393 

394 

395 
296 
39 ? 

398 

399 

400 
4 01 

402 
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404 
4 05 

406 

407 

408 

409 

410 

411 
4 12 
413 
4 14 

415 

416 
41? 

418 

419 

420 

421 

422 

423 
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425 

426 
42 ? 

428 

429 

430 

431 

432 

433 

434 

435 

436 
43 ? 

438 

439 


87 

IFF 


440 

00 

0 


43 

RCL 

DR51 

441 

00 

0 


51 

51 

+ DSPR 

442 

00 

0 


99 

PRT 

Print DSPR 

443 

00 

D 


87 

IFF 

If Flag 0 

444 

00 

0 


00 

00 

Is set, go 

445 

00 

0 


59 

I NT 

to 

446 

00 

0 


61 

GTD 

Subroutine 

44? 

00 

0 


49 

PRD 

INT; 

448 

00 

0 


76 

LBL 

otherwise 

449 

00 

0 


39 

CDS 

to to PRO 

450 

00 

0 


69 

DP 


451 

00 

0 


00 

00 


452 

00 

0 


02 

2 

U 

453 

00 

0 


03 

3 

n 

454 

00 

0 


03 

0 

o 

p 

455 
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0 


03 

3 

r 

456 

00 

0 


03 

3 

N 

45? 

00 

0 


01 

1 


458 

00 

0 


06 

6 

s 

459 

00 

0 


04 

4 


460 

00 

0 


00 

0 

o 

461 

00 

0 


01 

1 


462 

00 

0 


69 

DP 


463 

00 

0 


01 

01 


464 

00 

0 


03 

3 

Q 

465 

00 

0 


06 

6 

J 

466 

00 

0 


01 

i 

E 

467 

00 

0 


07 

7 


468 

00 

o 


03 

Z\ 

u* 

T 

469 

00 

0 


07 

r» 

( 


470 

00 

0 


00 

0 


471 

00 

0 


00 

0 


472 

00 

0 


69 

DP 


473 

00 

0 


02 

02 

Print 

HPN= 0 SET 

474 

00 

0 


69 

05 

DP 

05 

475 

476 

00 

00 

0 

0 


00 

0 


477 

00 

0 ! 


42 

51 

61 

59 

STD 
51 
GTD 
I NT 

0 + DR51 
Go to 

Subroutine 

INT 

4 78 
479 

00 

00 

0 

0 


76 

LBL 





89 

if 


005 

69 

DP 

List of 

53 

< 


099 

49 

PRD 

Subroutines 

24 

CE 

DS7R 

116 

44 

SUM 


65 

X 

* 

135 

60 

DEG 


73 

RC* 

DRtDROO > 

154 

70 

RFID 


00 

00 

169 

59 

INT 


54 

> 


343 

48 

EXC 


92 

RTN 


385 

87 

IFF 


00 

0 


395 

39 

CDS 


00 

0 


428 

89 
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00 
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00 

0 
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APPENDIX C 


THE SAMPLE PROBLEM AND ITS SOLUTION 
BY TI-59 PROGRAMMABLE CALCULATOR 


SEGMENT NO. li INPUT DATA AND COMPLETE DESIGN-POINT 

OUTPUT DATA . C-3 

SEGMENT NO. 2: OUTPUT DATA FOR SELECTED VALUES OF I DN 

AND FOR ANNUAL PERFORMANCE C-6 


C-l 


Sample Problem 


Statement, of Example 

Given (1) the required Input data describing a paraboloidal 
concentrator/Brayton solar thermal power plant operating without thermal 
or electric storage capability and (2) a histogram of the direct normal 
Insolation received at Barstow, California, In 1976, determine the design 
and annual performances of the power plant when connected to a conventional 
electric utility grid. 


PRESS OUTPUT/MODE 


ENTER 

i ' 

IDND 


1 . 

ENTER 

HOC 


0. 94 

ENTER 

HC 

0. 8S 

ENTER 

QRL 

0. 07445 

ENTER 

FSB 3 
0. 99 

ENTER 

ftLPH 


0. 92 

ENTER 

HP It 


0. 271 

ENTER 

HEGT 


0. 95 

J ! 


C-3 


COMMENT 


ENTER 


PRECEDING PAGE BLANK NOT FILMED 






ENTER 


PRESS 


OUTPUT/MODE 


COMMENT 


ENTER HPflR 

— 1 

0. 98 


ENTER flNPS 


0. 


-11. 484 

ao 

1. 


98. 721 

a l 

2. 


-303. 71 

a 2 

3. 


532. 7 

a 3 

4i 


-526. 26 

a 4 

5. 


275. 56 

a 5 

w » 

-53. 527 

a 6 

ENTER HOURS AT 

I UN 

0. 025 

-IDN 

251 • 

=HRS 

0. 075 

-IDN 

174. 

=HRS 

0. 125 

= IDN 

131. 

=HRS 

0. 175 

-IDN 

109. 

=HRS 

0. 225 

= I D N 

111. 

=HRS 

0. 275 

-ID N 

122. 

=HRS 

0. 325 

-IDN 

94. 

=HRS 

0. 375 

-IDN 

92. 

-HRS 

0. 425 

®IDN 

96. 

1 

=HRS 

1 _ 







ENTER 


PRESS 


OUTPUT/MODE 


COMMENT 


0. 475 
ill. 

=*IDN 

-HRS 

0. 525 
141. 

= IDH 
-HRS 

0, 575 
129. 

= IDN 
«HRS 

0. 625 
143. 

= I DH 
=HR8 

0. 675 
193. 

= IDH 
=HRS 

0, 725 

22ft, 

= I DN 
-HRS' 

0. 775 

O '“i © 

V i J . 

= IDN 
-'•IRS 

0. 825 
345. 

= IHH 
=HRS 

U » b { 5 
489. 

= IDN 
=HRS 

0. 925 
595. 

= IDH 
-HRS 

0. 975 
508. 

= i n h 

=HRS 

1. 025 
60. 

= ihn 
«HRS 


ENTER NPDL 



b. 

0. 94 

HDD 



. 8315978723 

-HRD 



0. £7 1 

HPD 



0. 931 

=HEP 



. 1972241963 

=HSD 



ENTER SECOND 

CARD 







ENTER 

PRESS 

OUTPUT/MODE 

COMMENT 



0. 025 

-I BN 

.000< I ON <.005 

— 

- 

0. 88 
. 9361702128 
0. 022 
. 0234042558 

- 

HC 

.HCN 

QC 

QCN 



-2, 473290909 
HRN=0 SET 

0. 

0. 

- 

HR 
HRN 
QR • 

QRN - 



“3. 112164 
HPH=0 SET 

0. 

0. 


HP 

HPN - — 

. QP - 

qpn 


- 

0. 

0. 

0. 

0. 

- 

. HS ... 
HSN 

qs 

qsN 



0. 075 

= IBH 

,005< IDN< .100 



0, 88 
. 9361702128 
U, 066 
0. 070212766 

-0. 217230303 
HRN=0 SET 

0. 


Note: Output 
parameters at 
each IDN appear 
in the same 
order as that 
shown above for 
IDN = 0.025 



U r 

-3. 112164 
HPN=0 SET 

0. 

0. 


Note: Because 

HEP = HEPD = 
a constant, 
qsN = qPN 
always. 



0. 

0. 

0. 

0. 


i 
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ORIGINAL PAGE IS 
OF. POOR QUALITY) 


ENTER 


OUTPUT/MOO E 


COMMENT 


0 . 125 


= I DN 


. 100 < IDN < .150 


. 9361702128 
0 . 1 1 

. 1170212766 

. 2339818182 
. £813641376 
0. 025738 
. 0329255906 


Note: This is the 
smallest tabulated 
value of IDN for 
which thermal 
output power 'from 
the receiver"1s 
positive -and, 
therefore, for 
which HRN >0. 


- 2. 360145651 
HPN-0 SET 

0. 

0. 


0 . 425 


= :dh 


.400< IDN 450 


0 . 88 
. 9361702128 
0 . 374 
. 3978723404 

. 7117358289 
. 8558653798 
0. 2661892 
. 3405251618 


Note: This is the 

smallest tabulated 
value of IDN for 
which net power 
is produced by the 
power conversion 
subsystem and, 
therefore, for 
which HPN>0. 


. 0921677418 
. 3401023681 
. 0245340574 
. 1158134139 


. 0537440176 
. 2725021504 
. 0228412075 
. 1158134139 


c- 







ENTER 


PRESS 


OUTPUT/MODE 


COMMENT 


1. 025 


Q oo 
■ oo 




o l I 


02128 


0. 902 
. 9595744 681 


. 828261 1973 
0. 995987684 
0. 7470916 
. 9557243041 


. 2708914533 
. 9995994587 
. 2023807293 
. 9553414971 


. 1838209356 
0. 932040485 
0. 188416459 
. 9553414971 


. 7862375972 


2489081067 


2847. 9 


U .16 0 3 3 3 8 4 3 


-I BN 


: HRfl 


= IE I 


=HSfl 


ORIGINAL PAG~ K 
OF POOR QUALITY, 


1.000 < ION <1.050 


c-a 





APPENDIX D 


USER INSTRUCTIONS FOR LOADING AND RUNNING 
THE TI-59 PROGRAMMABLE CALCULATOR PROGRAM 


» ■» Kb > L J & .. 


c\j co vr in m 

rtJ fO ro fTJ fO 



tn in in lo to ininintn intninin to in in in in 
cvji^cvjr^roi^cMr^rgr^cvjr^cvj oj n m n 
oor-r*(vjcvjfOfo<r<Mnnio <d r^^coco 


il II ii It ii n 

23 
ci s 


= o s 

tn 

u 

c Dc 
O 
2C 


S 8 
8 C 


II 11 II 


II II II II II II 


8 8 

c c 


8 K 
s a 


8 8 
8 8 


8 C 

8 8 


8 8 8 8 8 


c 

o 

ai 

n i 
c O 

•f- 

-u flj 
C O) 

a a 


“ s ?- ? ». a ?..? fi .s a 3 s « ,«-.s n a_ °_~ « » - 


P 4J D) 

V) c c 
•r- OJ »f- *u 
r- E 4J <U 

5i U) >> 

S dJ'r O ■ 

LO r— r~ 4-> 

D. C 

O. J E E* OJ 

°S2 s g, 

a) a> in qj 

a) O) Or™ 10 

c i. a) 

T- axi v) 
% m »r» 

♦ ♦ O r* r- JC 
0> r— 4-> 

4-> r— * L 

op op c 



(/> 

c 

o 


u 


3 


to 


C 

w 

0) 


w> 





OUTPUT/MODE { 


O 

z: 

o 

C£ 

UJ 

H 

2 : 

UJ 

X 

X 

X 

X 

ENTER HCD j 

>- 

X 

X 

X 

ENTER HC 

XXXZ 

ENTER QRL 

XXYX 

ENTER FSBS 

! AAXX 

x 

D. 

* 

cc 

UJ 

£ 

LU 

ZAXX 

ENTER HPD 

1 XZ XX 

ENTER ANPS ! 

~ 

» 

o 


PRESS 

















— 

1 


i 

Irst 

to 

a; 

to 

cc 


s/a 


R IS 


R/S 


CO 

cc 


to 

cc 


to 

cc 






















cc 

LU 

1- 

Z 

ID 



mi 


XXXY 


XXXZ 


XXYX 


XX YY 


XX YZ 


XXZX 


1 


! PROCEDURE 

Start program segment no. 1 


Enter IDND 


Enter f!CD 


Enter HC 


Enter QRL 

j 

Enter FSBS 

i 

| 

1 1 

! i 

i 

i 

Enter ALPH j 


Enter HPD 


— 



STEP 

H 


OJ 


ro 




LO 

L 

- 




CO 





D-3 


PFJCCEDING PAGE BLANK NOT FILMED 





hours continues at each 0.05 
increment of IDN to IDN = 1.025 
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7 is stored in 03 a 
»step 177 of seg. 
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